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I.  INTRODUCTION 


The  signals  at  the  input  to  the  difference  and  sum  channels 
of  a monopulse  receiver  are  (in  complex  envelope  notation) 

D(t)  - [gp  (0)  gs  (0)  A + Nd]  eiut  (1) 


S(t)  - [gs  (0)  g s(9)  A + Ng]  eltjt  (2) 

Where  gp(0)  and  gg(0)  are  the  one  way  voltage  gains  (for  off-boresight 

angle,  9)  of  the  difference  and  sum  antennas;  A is  the  voltage  return  at 
an  isotropic  antenna  (it  includes  transmitter  power,  target  range,  target 
backscatter  coefficient,  etc.);  and  N^(t)  and  N<,(t)  are  zero-mean  Gaus- 
sian receiver  noises.  It  is  assumed  in  Equations  (1)  and  (2)  that  the 
target  has  been  illuminated  by  the  sum  pattern  and  received  by  both 
difference  and  sum  patterns. 

Alternatively,  the  following  can  be  written: 

Dr  ■ Real  D(t)  - x cos  u>t  - y sin  ut  (3) 

SR  ■ Real  S(t)  « u cos  wt  - v sin  u>t  (4) 

where  x,  y and  u,  v are  the  in-phase,  quardrature  components  in  the 
difference  and  sum  channels,  respectively. 


The  monopulse  ratio  is  defined  as 

/ W* 

xu  + yv  LPF 

2 . 2 ■ f_ 2.„ 

u + v I S_dt 

ih  R 

Clearly,  when  there  is  a single  target  and  no  noise, 

gn(e) 
ri  * H<fi) 


A „ , 

- Real 


\S(t)f 


(5) 


(6) 


Since  the  right  side  is  a monotonic  function  of  9 (to  avoid  ambiguities) 
in  the  region  of  interest,  measuring  r^  yields  9 , the  angle  of  the  tar- 
get with  respect  to  boresight.  This  result  is  independent  of  target 
amplitude  fluctuations  which  are  manifest  in  A. 


When  noise  is  present,  the  following  Is  obtained: 


(x  + n ) (u  + n ) + (y  + tiy)  (v  + n ) 

1 — 2 2 

(u  + nu)  + ( v + n^ 


(7) 


where  the  over  bar  represents  a statistical  mean  caused  by  the  presence 
of  the  target.  Because  the  phase  reference  is  arbitrary,  the  target 
phase  can  be  used  as  a reference  i.e.,  the  target  is  set  completely  in 
the  in-phase  components;  then. 


(x+n)(u+n)  + n n 
x x u y v 

— 2 2 
(u  + n ) + n 

u v 


II.  APPROXIMATE  ANALYSIS 

At  this  point  an  approximate  analysis  is  usually  made  which 
goes  as  follows:  Since  the  signal-to-noise  ratio  is  large,  the  quadra- 
ture terms  (which  have  no  signal)  compared  to  the  in-phase  terms  can 
be  neglected: 


x + n 


r,  « 


1 u + n 


u 


(9) 


so  that  r^  is  the  ratio  of  two  Gausssians.  Again,  using  the  fact  that 
the  signal-to-noise  ratio  is  large,  i.e.. 


-2 

u 


» n 


(10) 
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Che  noise  in  the  denominator  of  Equation  (9)  compared  to  the  signal 
there  can  be  neglected;  thus, 


(ID 


and  because  n is  zero  mean 
x 


(12) 


l.e.,  the  estimate  is  unbiased,  and 


(13) 


—2  2 
Since  u is  the  total  sum  signal  power  and  n^  is  half  the  noise  dif- 
ference power,  [Equation  (8)],  the  accuracy  of  the  monopulse  measurement 
is  taken  as 


1 

°ri  ’ yrm 


(i4) 


This  approximate  analysis  linearizes  the  monopulse  ratio  and  states 
that  the  average  of  a quotient  equals  the  quotient  of  the  averages 
(of  numerator  and  denominator) . Clearly,  Equation  (12)  cannot  be  valid 
because  by  decreasing  nu  in  Equation  (9)  by  a small  amount,  e,  r^  is 

increased  more  than  it  is  decreased  by  increasing  by  an  equivalent 

amount.  Thus,  a bias  must  exist.  A further  critique  of  the  passage 
from  Equation  (9)  to  Equation  (11)  is  in  order.  The  validity  of  Equa- 
tion (10)  does  not  allow  the  assumption  that  nu,  the  instantaneous  value 

of  the  noise,  is  always  less  than  u.  After  all,  the  noise  is  Gaussian 
and  can  assume  all  values. 


In  the  exact  analysis  given  in  the  following  equations,  tha  bias 
will  be  made  manifest  and  no  assumptions  concerning  large  signal-to- 
noiae  ratio  will  be  Invoked. 
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III.  EXACT  ANALYSIS 


Equations  (1)  and  (2)  can  be  rewritten  for  the  situation  of 
multiple  targets  as 


D(t) 


W *s(e 


) A + N_  ] 
o'  n D I 


iu)t 


(15) 


N 


S(t) 


. n=*l 


gc(0„)  S(e)  A + N.  | e 
b n s n n S 


iujt 


The  noise-free  monopulse  ratio  is  defined: 
N 


a A 
r - Re 


r B 

n“l  n n 
N 

y b 

L n-1  " 


where 


A 2 

B “go  (0  ) A 
n °S  ' n'  n 


th 


is  the  sum  channel  return  from  the  n target  and 


A Sp^n^ 

^ 8S«V 


th 


(16) 


is  the  monopulse  ratio  associated  with  the  n target.  The  random 
phases  of  the  targets  are  included  in  the  complex  numbers,  A^,  and 

give  rise  to  glint,  but  in  this  analysis  they  are  treated  as  constants. 
As  distinct  from  the  single  target  case,  the  following  no  longer  applies: 


x . y 

™ m • 

U V 

The  following  auxilliary  variable  is  introduced: 

, . A . .2,2. 

t (r)  ■ xu  + yv  - r(u  + v ) ; 

»x,u,y,v 


(21) 
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\ 

/ \ 

\ 


/ 

/ 

/ 


then 


Pr.  (rx  £ r)  = Pr.  0)  , (22) 

thus  the  probability  density  function  (PDF)  of  r.  is  given  in  terms  of 
the  PDF  of  ^ by  1 

o 

Pr  (r)  ■ ^ j p (a)da  . (23) 

^ -OO 


Note  that  £ does 
The  Fourier 


not  involve  a ratio, 
transform  of  p^(a)  ie  defined  by 


00 

Fj,  (a j)  - J p^(a)  e_l0!uJda 

-00 


(24) 


and  the  Fourier  integral  theorem  states 


00 

P^a)  m Tn  f F^(w)e+la“d“ 

-00 


(25) 


It  appears  at  first  glance  that  to  find  p^,  is  needed  and  con* 

versely  to  find  F^,  p^  is  needed  — so  what  is  the  use  of  introducing  the 

Fourier  transform?  However,  Equation  (24)  Btates  that  F^(u)  is  the 

expectation  of  e This  may  be  calculated  over  the  domain  of  the 

random  variables  x,  u,  y,  v.  Thus, 


F £ (a) 


00 

M- 


2 2 

-iu[  xu  + yv  - r(u  + v )J  , . . . 

e Px,u,y,v  («,u,y,v)d*...dv. 

(26) 


The  in-phase  and  quadrature  noise  components  in  each  channel  are 
independent  and  identically  distributed.  Independence  between  the  noises 
in  the  difference  and  sum  channels  of  the  receiver  is  assumed  so  that 
the  joint  density  in  Equation  (26)  becomes  the  product  of  the  four 
densities 


Px(x) 


(27) 
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py(y) 


—i. e -to h2 

c 2 
In  a 2a 


- > 1 -(u-u)^ 

p»(u)‘  ^ 

Pv(v)  “ ~ZH  e ' 

S/Jnrb  2b 


?x,u,y,v(x’U»y,v)  * px(x)  pu(u)  Py(y)  Pv(v) 


Hence 


/‘-1 


-iu(xu-n»  ) 


Px(x>  Pu(u)  dx  du 


“toj(yv-rv  ) , . . . . 

e py(y)  Pv(v)  dy  dv 


The  double  integral  in  x,  u is  evaluated  by  completing  the  squares 
in  the  exponent  and  yields 


iv — + ^ ii-  r~  *i2  * i 


( )dx  du  ' 


. 2b2  IV*  u / I 2bz  /m b 

• • J>  b “2  - 2llr.hu  + 1 


^ • b w • 2ljr«bu  ♦ 1 


The  double  integral  in  y,  v is  obtained  from  Equation  (33)  by  replacing 
x,  u by  y,  v,  respectively.  Thus  j 


- £u2*2  + ^»2  + a2  ♦ v^j 


*,<“>  - £ s-s- 5 ; 

5 tVu  - 2i»ibu‘  + l 


,,  . U-2^  . “ 2^2  \ —2'.  — 2~  ”1  T —2^2  —2^2  —2 

..  2UbM  LVU  + v *2  )•  * tt  »1  * v '2J-Lu  «i  ♦ v '2  + ? ♦ 

2b2  r»Vu2  - 2K.hu  + ll 
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where 


s(r)  £ 

a 


A A b x 

si=i: 


A A b v 

8-  « — * 

2 3 V 


/ 


/ 

When  only  a single  target  is  present,  %-,m  s , . 

/ - ' 


/ 


/ 


In  the  absence  of  noise  tL-j  mraopuioe  ratio  i >'  given  by 

r m x u + v 1 
-2-2 
u + v 


(35) 

(36) 

(37) 


(38) 


Define  also 


^ A br 

a 


—2*  — 

* U + V Sj 


u2  + 


and 


(39) 


~2~2  -tf'2  -2-2  - - 2 

-2  A u + v s2  ^ uV(8l  - s2)Z 


u2+  v2 


- %*■  +■ 


(a2  + v2)2 


& 0 


For  a single  target  t - s. 


(AO) 


•lmpllflei'co  aW  E<IU*t,0n'  <38)-  <39>-  “4  <«».  (34) 


F?(a>)  - 


■X(t2  + l)e-«^la^  S+B)  - (t2  + 1) 
a b u - 2lsabu,  + l 


2.  2 2 

a U {Jj 


J jL 


(41) 


f • . 


where  the  sum  channel  signal-to-noise  ratio  has  been  introduced- 


x - °2  ± v2 
2b2 


(42) 
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Since  this  is  not  a function  of  s - s unless  s ■ 0,  or  from  Equation 
(39)  'X/TT  - 0 ■ y"/v,  it  was  concluded  that  unless  there  is  a single  target 
on  boresight,  there  will  be  a bias  in  the  estimate  of  r^.  Substituting 

Equation  (25)  into  Equation  (23),  the  following  is  obtained: 


Pr,  <r>  7 


■k  f T.  V“>  0" 


dt; 


(A3) 


or 


Pr  (r) 
rl 


o oo  la  r 

b f i f d . , \ ab  Sl(j 

" J = J “Fc  U)e  s 


d£ 


(44) 


Since  F^  (u/ab)  is  given  by  Equation  (41),  the  differentiation  performed 
and  the  following  is  obtained: 


P_  (r) 
rl 


2b  -X(t  +1) 

a e 


-oo  _00  L <jj  -2isw  + 1 


2 

.x  2 MM  -t,.?) .._Ct_2  +J± 

2 

u)  - 2iscj  + 1 

ft 

(m2  - 2isu>  + l)2 


(45) 


! 


In  carrying  out  the  u integration,  the  factor  iio/ab  corresponds  to 
differentiation  with  respect  to  £;  this  combines  with  the  subsequent 
integration  over  £ from  ■»  to  0 so  that  the  result  of  the  to  integration 
needs  to  be  evaluated  only  at  £ ■ 0.  Thus  in  terms  of 


G(t) 


00 

A 1 ft 

" T«  J 


( 2lu>(t28  + 8)  - (t2  + 1)\ 

\ U2  - 2isu)  +1  / 

— — du 


(uj2  - 2isu  + l)2 


the  PDF  of  the  monopulse  ratio  becomes 


2b  a -x'?2+1>  r<l-X*2)G(T)-K;’<T)l 

»r  <*>  s — L Jt«1  . 


(A6) 


(47) 
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I 


V 


n *.1*  CTdfelut0  evaluate  G(*>.  the  integrand  is  written  as  a product 
n ter.i.s  of  the  ssaential  singularities  at  ™ 


v , 2 - if  a ±Js2  + 1) 


where 


oo  . 1BT 

, J_  f e “'“I  e u_w2 

2n  J ( \2  * 7 72  a° 

-o=  vu"wl/  (<*>“u2) 


{ [4  +J^)  * i]2  * <?2  - i2>  (.  +/!Ti)2| 


2 7X7 

vs  +1 


2 fl 

V B +]. 


{[;(,  -yPTi),  .]*  + (?  l (.  . 


tinf:e  W1  liAa  in  ’;he  u^'er  h91"  P^ne  and  ^ lies  in  the  lower  half 

by  “"Pletln*  the  • 


i-  / -g  — - M1  „iu>t  j 1 + aent 
2 t f , TI  e dw  

:C.  ' “J  2 


t-D  .^)‘\(^6nr) 


00  -J£J. 

/•  i.i  - <■ 


_1_  | e u " w 

2n  7 "7 — 

(U  - u, 


p^7 


where  ^ denotes  the  modified  Bessel  function. 


Thus  using  Farseval's  theorem,  Equation  (49)  becomes 


G(t)  =*  / e 


-2  /m  t (^1/2  ^1/2 1[(2  ^ ) ,t  (2 


Expanding  the  Bessel  functions  into  power  series  and  employing  the 
Cauchy  product,  the  following  is  obtained: 


00  rr—  00  T n 

■/  ‘2  1 

0 n«0  . m-0 


(n  - m) ! ( 1 + n - m)  I 


m!  (1  + m) ! 


2+n  . 

t dt 


Termwise,  integration  yields 


G(t)  - 


t=-  I <* 

9 9 


+ n)! 


(2  /s  +1)  n-0 


1* 

2 ml (1  + 


m)  I (n  - m)  I (1  + n - m)  I 


for; )‘ ' 


Since  p/a  is  independent  of  X,  this  yields  G(t)  and  hence  G'  (*r)  and 
hence  the  PDF  of  the  monopulse  ratio  as  a power  series  in  a which  is 
proportional  to  the  signal-to-noise  ratio.  This  is  most  inconvenient 
for  large  x.  To  obtain  a representation  in  closed  form,  it  is  noted 
that  the  polynomial  in  parentheses  is  a degenerate  hypergeometric 

function*,  namely: 


Y 

/.  m!  (1  + m)l  (n  - ra)!  (1  + 


n - m)J 


n!  (1  + n)! 


2* l(-n,  -n  -1;  2;  p/a) 


Abramowitz,  M.  and  Stegun,  I.  A.,  Handbook  of  Mathematical  Functions. 
National  Bureau  of  Standards  Appl.  Math.  Series  Vol.  55,  March  1965 
(Formula  15.1.1). 
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This  is  expressible  in  terms  of  a Legendre  function  of  the  first  kind; 

2 

the  following  quadratic  transformation  is  employed: 

2Fl(-n,  -n  -1;  2 ; z)  = z‘1/2  (1  - z)1+n  P_l  (58) 

1+n  U ' z/ 

where  Pj^  is  the  Legendre  function  of  the  first  kind  of  order  -1  and 

degree  1 + n.  This  is  valid  whenever  |a*g(l  - 2) ) < it  and  z does  not 
lie  between  0 and  Since  0/a  is  positive.  Equation  (58)  holds  when- 
ever z (*  0/a)  < 1.  Assuming  that  a > 0,  then  Equation  (58)  holds. 
Using  the  relation  between  the  Legendre  functions  of  the  first  kind  wit 
negative  order  and  the  Legendre  functions  of  the  first  and  second 

3 

kind  with  positive  order,  the  following  is  obtained: 


p-1  (.L±-£) Si p1 

rl+n  \1  - zJ  (2  + n)J  rl+n  \1  - z/ 


(59) 


thus 


SSM 


m 


m"0 


m!  (1  + m)!  (n  - m)J  (1  + n - m)| 


1 l pi  (a  + 9 \ 

(2  + n)t  (1  + n)l  rl+n  \a  - 0/ 


(60) 


Since  for  w > 1,  the  following  is  obtained: 
2 


,1 

' 1+n 


p:._  <w)  - (w‘  - 1)1/2  (w) 


(61) 


where  (w)  is  the  first  derivative  of  the  Legendre  polynominal  of 

degree  l + n,  G(t)  becomes 


G(T) ' 2 jrhyi  (^i)(^i^=) 


n*0 


ibid.  (Formula  15,4.14). 
*Ibid.  (Formula  8,2.5). 


(62) 


Ibid,  (Formula  8.6,6). 


Since  G(t)  is  invariant  with  respect  to  an  interchange  of  a and  6 
[Equation  (54) ] , if  p/a  $ 1,  then  a and  p can  be  interchanged.  Thus 
a-p  in  Equation  (62)  may  be  replaced  by  ja-p|  and  Equation  (62)  will 
hold  for  all  a and  p.  Returning  to  Equation  (47),  the  PDF  of  r,  is 
written  as 


/ \ be 
\ (r>  ■ 2a 


A 2 

-x(t  +1) 


m=l 


(m  'Xt  ) „ (1)  m-1 


m! 


? ' ' t 

m(z) 


(63) 


where 


£ lls  , 

|a  - p| 

£ JoL-^fe] 


(64) 

(65) 


Since  PQ  (z)  is  a polynominal  of  order  zero  inz,  the  lower  summation 
limit  in  Equation  (63)  is  extended  to  zero. 

A generating  function  for  the  Legendre  polynomials  is  given  by5 

I V*>£-*tZlo  (66> 

m=0 


hence  Equation  (63)  becomes 

‘2 


pri  (r)  " 2a 


e “X(t  +1)  1 J_l  J_  v^lTtz  ( /2  _ Y) 

(s2+d  3/2  c dzr dt  n °v  / 


Carrying  out  the  t differentiation,  the  following  is  obtained: 


b e i d i *2 

Vr)'S  (.2  + 1)W2  t * Y 


o(^) 


+e 


tz 


{^)h  (*/^ i) 


+ 

(68) 


JIbid.  (Formula  22.9.13). 
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where  the  following  relation^  has  been  used: 


3c  lo  “ h <■=> 


Now  let 


..  A . a + 

W = Cz  = ~ - 


2 ./S  + 1 


so  that 


. -X(t  +1)  ( w.  ~2.  / 

Pr  <r>  a s vT  £ £ (W'Xt  } M* 

rl  2a  (s2  + 1)3/2  dw(  V 

+ ” jj-  - t2  i^Tw^Tt2)!  . 

The  following  function  is  introduced; 

H(m,k,w)  = eW  wm  ^v42  - t2)  1^  (^/w2  - t2 


(w-xt2)  IQ  f / w2  - t2  ) + 


/ > be  ~X^£  +1^  d ( \ 

pri<r)  “ 2a  TT^.y/i?  35  }H<l»°»w>  ' ^tTKO.O.w)  ^H(0,l,w)| 


Since 


[x\(x)]a  *k*llk-1(x)  ’ ( 

differentiation  of  the  composite  function  H(m,k,v)  yields 
d 

H(tn,k,w)  =»  H(ra,k,w)  + mH(m  - l,k,w)  + H(m  + 1,  k - l,w)  . 


Ibid.  (Formula  9.6.27). 
2Ibid.  (Formula  9.6.28). 
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Employing 


8 


!k(*)  ■ 


(76) 


and  collecting  terms,  the  following  Is  finally  obtained: 


. -x(t2+l)+w 

Pr  * 2a  2 3/2 

rl  Za  (sZ  + 1 y/l 


+ I*L 


2 2 
wxt  - t 


JO 


[(2w+l-Xt2)  IQ  U2  - t2)  + 

i h(JT77)]  . 


(77) 


This  Is  closed  form  solution  valid  for  any  monopulse  radar,  and  It  Is 
suitable  for  computation.  It  yields  the  PDF  of  the  monopulse  ratio 
(which  as  pointed  out  In  the  text  is  biased) , as  a function  of  the  mono- 
pulse ratio,  r.  If  p (0)  Is  desired,  the  monopulse  ratio  as  a function 

®1 

of  the  off-boreslght  angle  6 , one  merely  starts  with  the  one-way  voltage 
antenna  difference  and  sum  patterns  and  calculates 


r(Q) 


* ^>(9) 


8g(9) 

(Figure  1)  and  Its  derivative,  then 

■>„  <«)  - Pr[r«)l  |f  (0)|  . 


(78) 


(79) 


There  are  many  papers  in  the  literature  which  attempt  to  analyze 
the  performance  of  a monopulse  radar.  All  but  one  confine  their  atten- 
tion to  a linear  r(0)  characteristic  and  assume  large  slgnal-to-noise 
ratio..  None  come  to  grips  with  the  nonlinearity  inherent  in  the  mono- 
pulse ratio  and  hence  do  not  arrive  at  a bias  or  a PDF.  They  are  limited 
to  approximating  the  variance  of  the  PDF.  In  order  for  the  sceond 
moment  to  exist,  i.e.. 


(r)  dr<  « 


(80) 


8.  Ibid.  (Formula  9.6.6). 
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the  PDF  must  fall  off  as  |r|  «,  faster  than  r . For  large  |s|, 

a and  £ approach  a constant  times  s,  and  w and  t approach  constants. 

_3 

Thus  p (r)  falls  off  only  as  fast  as  r . It  is  concluded  that  the 
rl 

second  moment  does  not  exist  and  hence  that  the  variance  is  infinite  for 
any  signal-to-noise  ratio.  However,  the  distribution  has  a spread  which 
may  be  calculated  by  comparison  with  that  of  a Gaussian  PDF  as 


(81) 


where  R^  and  R£  are  given  by  the  solution  of  the  equations 


(82) 


(83) 


Since  p (r)  given  by  Equation  (77)  is  convenient  to  program,  R and  R, 

rl  l . 

are  readily  calculable. 


The  only  analysis  which  does  not  ignore  the  fundmental  nonlinearity 

9 

in  the  monopulse  ratio  is  one  in  the  Soviet  literature  . The  authors 
employ  a brute  force  transformation  of  variables,  and,  after  integrating 
over  three  of  the  four  variables  involved,  succeed  in  obtaining  the 
following  expression  for  p (r): 

rl 


(r+ga)  „ _b_ 
*r^  b 2a 


(-D*(2  -V 


(2k  + 1)1  j 
(2k)! 


9 

Alexandrov,  V.  G.  and  Fedosv,  V.  P.,  "Statistical  Characteristics  of 
Signals  at  the  Output  of  a Monopulse  Radar  Receiver  with  Non  Fluctuating 
Input  Signals,"  Academy  of  Sciences.  U.S.S.R. . Radioelectronics.  Vol.  17, 
No.  4,  April  1974  (in  Russian) . 
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and  6^q  is  a Kroneker  delta,  1^  la  the  modified  Bessel  function  of  order 

k and  jF^  is  the  confluent  hypergeometric  function.  Since  each  term  in 

the  infinite  series  involves  a different  Bessel  function  and  hypergeo> 
metric  function,  their  expression  is,  to  say  the  least,  unwieldy. 

The  procedure  to  be  followed  to  calculate  the  PDF  of  r^  is  to  choose 

the  target  location  or  locations  arbitrarily  and  assign  values  to  x,  y, 
u,  v,  a,  b.  Then  Equation  (42)  gives  X.  For  a selected  value  of  r, 
s(r)  is  calculated  from  Equation  (35),  fc.  and  s_  from  Equations  (36) 

and  (37).  Then  s and  t are  uvallable  from  Equations  (39)  and  (40). 
Equations  (50)  and  (51)  give  a and  (3  while  Equations  (65)  and  (70)  give 
t and  w.  These  values  are  employed  in  Equation  (77)  and  as  r is  varied, 
the  PDF  of  r^  is  computed. 
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